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The spontaneous formation of vortices is a hallmark of collective cellular activity. Here, we 
study the onset and persistence of coherent angular motion (CAMo) as a function of the number 
of cells N confined in circular micropatterns. We find that the persistence of CAMo increases 
with N but exhibits a pronounced discontinuity accompanied by a geometric rearrangement of cells 
to a configuration containing a central cell. Computer simulations based on a generalized Potts 
model reproduce the emergence of vortex states and show in agreement with experiment that their 
stability depends on the interplay of spatial arrangement and internal polarization of neighboring 
cells. Hence, the distinct migrational states in finite size ensembles reveal significant insight into 
the local interaction rules guiding collective migration. 


PACS numbers: 87.18. Gh, 87.17.Jj,87.18.Fx,87.17.Aa 

The ability of cells to coordinate their motion is es¬ 
sential in various biological contexts, notably morpho¬ 
genesis m and tissue repair mis]. In recent studies, 
monolayers of Madin-Darby canine kidney (MDCK) cells 
have been investigated as model systems for collective 
behavior in living systems. Remarkably large scaled cor¬ 
relations and swirls in cell migration have been observed 
and characterized using image correlation and traction 
force microscopy techniques [SHS]. These emergent pat¬ 
terns and correlations were attributed to cell-cell cou¬ 
pling, and mechano-transduction mediated by the force¬ 
generating cytoskeleton. In fact, dynamic self-ordering 
into streaming patterns and vortex states appears to be 
rather generic in assemblies of (self-)propelled objects. 
They are well known in active systems as diverse as 
driven biopolymers in motility assays [SHII], bacterial 
colonies [I2HI11, and driven granular media [T5HT7] . Su¬ 
perficially, these phenomena may be attributed to a ten¬ 
dency of neighboring objects to align their direction of 
motion, as suggested by flocking models [18]. However, 
upon closer inspection, there are many important quali¬ 
tative differences between all these systems and to date 
quantitative theoretical models are largely lacking. 

For cell assemblies, the challenge is that mechani¬ 
cal and biochemical interactions between cells as well 
as internal organization of cells are complex m. and 
therefore parameter control is limited. Recent progress 
in understanding collective behavior of cell assemblies 
has been fueled by micropatterning techniques which 
enabled well-controlled in-vitro experimental systems. 
These techniques have been used to study static adher¬ 
ence and intracellular cytoskeleton organization of indi¬ 
vidual cells in defined geometries [3D]. Importantly, geo¬ 
metrical confinement of cells into micropatterned circles 
has been found to induce persistent rotational motion 
for systems ranging from two cells m to large assem¬ 



FIG. 1. (a) Array of MDCK cells seeded on circular mi¬ 

cropatterns. (b) Circular patterns occupied by 2-8 cells. Cir¬ 
cle size increases in such a way that the average area per cell 
is constant at approximately 830 pm^. Nuclei are labeled in 
blue, (c) Schematic of 4 cells rotating within a circular field. 

blies im [23] (23] . There is general consensus that on a 
macroscopic scale collective cell migration is to a large 
degree generic and can be explained by different classes 
of theoretical models including flocking models UHlIlil- 
HZ], cellular Potts models [25H52] . and phase field mod¬ 
els [331 [31] ■ However, the mechanisms underlying the 
emergence of vortex states are still poorly understood 
and its relationship to single-cell properties remains un¬ 
clear. In particular, a systematic study of the emergence 
and stability of small-scale vortex states and the dynamic 
disorder-order transition leading to the emergence of col¬ 
lective migration as a function of the number of involved 
cells has not been carried out so far. 

Here, we investigated the emergence of collective rota¬ 
tional motion in small circular micropatterns as a func¬ 
tion of the number of cells (Fig. 1). The physical sys¬ 
tem consists of arrays of circular helds containing 2-8 
cells. Cell density is kept constant by increasing field 
size in line with cell number. We found distinct transi¬ 
tions between states of disordered motion (DisMo) and 
states of coherent angular motion (CAMo). Further- 
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more, the survival time of the coherent state tends to 
increase with increasing cell number, but shows a pro¬ 
nounced drop between 4 and 5 cells, where the geometric 
cell arrangement changes from a conformation without a 
cell in the system center to one including a centered cell. 
Employing a computational model, based on the cellular 
Potts model (CPM) [3S1121], which we extended to in¬ 
corporate internal polarization and cell-to-cell mechano- 
transduction [3S] , we reproduced and explained these fea¬ 
tures. Thus, the experimentally observed gradual transi¬ 
tion with increasing system size from predominantly er¬ 
ratic motion of small cell groups to directionally persis¬ 
tent migration of larger assemblies is captured by the 
theory, underlining the role of internal cell polarity in 
the emergence of collective behavior. 

Micropatterns of the extracellular matrix protein fi- 
bronectin separated by PEGylated cell-repelling areas 
were fabricated using a plasma-induced patterning ap¬ 
proach. Parts of a culture dish (Ibidi) were covered with 
a polydimethylsiloxan (PDMS) template of the desired 
pattern. Exposed parts were treated with 02"Plasma in a 
plasma cleaner (electronic diener) and overlaid for 30 min 
with 1 mg/mL PLL(20)-g[3.5]-PEG(2) (SuSoS). After¬ 
wards, the template was removed and the whole surface 
was briefly exposed to a 50 ^g/mL solution of fibronectin 
(Yo Proteins). MDGK cells were seeded on the struc¬ 
tured surface and placed in a temperature-controlled en¬ 
vironmental chamber on the microscope stage. Arrays of 
circles were designed with increasing sizes to accommo¬ 
date 2-8 cells (Fig0. To ensure constant cell density of 
830 ^m^/cell, for each pattern size, only fields containing 
the appropriate number of cells were selected for analy¬ 
sis. Nuclei were stained using Leibowitz L-15 medium 
(c-c-pro) containing 15 ng/mL Hoechst 33342 (Invitro- 
gen). Time lapse movies were recorded at a rate of 6 
frames/hour over 50 h using an iMIC automated micro¬ 
scope (TILL Photonics). Individual nuclei were tracked 
using in-house image analysis software. 

A typical array of circular adhesion sites occupied 
by MDGK cells is shown in Fig. [^a). Cells exhibit 
spontaneous collective rotation within the circular areas 
(Fig-H c)). Periods of CAMo are seen to be interrupted 
by intervals of DisMo, after which rotation in an arbitrary 
direction is resumed (for movies see [36jl. Increasing the 
system size cell by cell (Fig. [^b)), we studied collective 
rotation as a function of cell number. For each cell i, 
the center of the nucleus was tracked and recorded in 
polar coordinates, and the individual angular positions 
(^i(t) were calculated (Fig. [^a)) (for a detailed descrip¬ 
tion see section S2 [35]). Typical time courses of (pi{t) 
for a system of 7 cells are shown in Fig. [^b). To filter 
out small fluctuations which result, for example, from 
displacements of the nucleus with respect to the geomet¬ 
ric center of the cell, we calculated the system angular 
velocity ^lN(t) as the mean over the individual angular 
velocities of the A^-cell system smoothed over a number 



FIG. 2. (color), (a) False-color fluorescence image of the nu¬ 
clei of seven cells within a circular micropattern. For each 
nucleus i, the angular position was evaluated with re¬ 

spect to the circle center, (b) Angular positions of each 
cell (in colors corresponding to the nuclei in (a)) and normal¬ 
ized total angular velocity The classification threshold 

of = 1/4 is indicated by the red dashed line. Periods of 
DisMo are highlighted by gray shaded areas, (c) Probability 
distribution of the mean angular velocity |fljv| for systems 
containing 2 to 8 cells. The distributions are fitted by a sin¬ 
gle Gaussian (green) and a mixture of two Gaussians (dashed 
red). The deviation between the two curves reveals a local 
maximum at |f2jv| =0. (d) Log-log plot of the angular MSD 
of CAMo (green) and DisMo (red) and its error for assemblies 
consisting of 8 cells. For other cell numbers see Fig. S4 |36 |. 

of frames u/ taken in discrete intervals of = 10 min: 

N t+nf 

^N{t) = ^ 2 ~ ~ 

'v 7 i=l T-t 

We chose n/ = 9 as the best trade-off between smoothing 
of fluctuations and temporal resolution [3H]. For all N, 
the probability distribution P(nAr) displays symmetry 
breaking into clockwise and counterclockwise rotations. 
Both directionalities are almost equally represented, with 
a small bias towards clockwise rotation (see Fig. S2 [36)1. 
Similar chiralities have been reported before miMiiin]. 

To distinguish periods of GAMo from periods of 
DisMo, we analyzed the probability distribution P(|r2Ar|). 
It was found to be approximately Gaussian (Fig. Sc)). 
The maximum, fijv, as well as the standard deviation, 
CTjv, decreased with increasing cell number, displaying 
an almost constant coefficient of variation crjv/^Af = 
0.74 ± 0.13. At Dtv = 0, P(|D 7 v|) exhibits a weak 
second maximum, indicating a state of disordered, i.e. 
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FIG. 3. (color online). (a) Survival function 5']v(t) = 
Pn{T > t) of CAMo and DisMo states. Insets show cor¬ 
responding log-lin plots. Exponential fits are indicated by 
dashed lines (for other cell numbers see Fig. S5 [36]). (b) 
Persistence time r as a function of cell number, derived from 
experiment and theory. Error bars indicate confidence bounds 
of 99% within the fits, (c) Peak positions Djv of the distri¬ 
bution of the angular velocity P(fljv) from experimental data 
and theory. Error bars indicate the standard deviation. 


non-rotating, motion. Introducing a normalized variable 
^N{t) ■= |nAr(t)/riAr|, we defined a common threshold 
for all fV at = 1/4, so that for ^Ar(t) < a migration 
state is classified as DisMo and for ^Ar(t) > as CAMo 
respectively (Fig. S2 [36]). (As discussed in section S3 
[36j , an alternative approach to identify collective motion 
gave the same results). To verify that these two states are 
distinct in their migrational behavior we calculated the 
angular mean squared displacement (MSD) during each 
state, MSD(t) = ([(i^(t))Ar-(</5(0))Ar]^)states, where t = 0 
signifies the starting point of an interval. Averages were 
taken over all N cells within a given system as well as 
over all observed intervals of CAMo or DisMo, denoted 
by (.. .)jv and (.. .)states, respectively. Consistently, the 
MSD of CAMo shows a slope 2 in a log-log plot, indicat¬ 
ing ballistic angular motion for all cell numbers, while the 
MSD of DisMo exhibited diffusive behavior (Fig.j^d)). 

Next we evaluated the lifetimes of the CAMo and 
DisMo states. Fig. [^a) shows the survival probability 
SN^t) = Pn{T > t), i.e. the fraction of CAMo/DisMo 
time periods T exceeding t, based on a sample size of 
over 600 systems (see Table S2 [36]). We found that 
the survival probabilities of both states decay exponen¬ 
tially, cx suggesting that the stochastic pro¬ 

cess underlying the emergence and collapse of both states 
is Poissonian. The persistence time t of the coherent 


state increases with increasing cell number, but exhibits 
a pronounced discontinuity between systems containing 
4 and 5 cells (Fig. |^b)). 

To further explore the mechanism underlying the dis¬ 
continuity in persistence time, we monitored the spatial 
arrangement of cells within the pattern. Fig. |^a) shows 
the relative positions of cells with respect to a reference 
cell. In systems containing up to 4 cells, the cells are pre¬ 
dominantly arranged in topologically equivalent positions 
in the outer regions of the circle. In this configuration, 
cells in the state of CAMo follow each other in a closed 
circle. As the number of cells increases to 5, the packing 
geometry changes abruptly to a conformation in which 
a single cell is located at the system center. To connect 
this topological transition to the observed decrease in the 
persistence of the CAMo state, intrinsic cell properties 
have to be accounted for. It is generally assumed that a 
migrating cell is highly polarized with respect to protein 
distribution and cytoskeletal organization [T9l[4T]. In ad¬ 
dition, since neighboring cells are coupled mechanically 
by cell-cell adhesion, a cell obtains directional guidance 
cues from adjacent cells. This coupling suggests that ad¬ 
jacent cells tend to align their direction of internal front- 
rear polarization. Hence, a ring-like arrangement, as seen 
for 2-, 3-, and 4-cell systems, naturally provides a stable 
conformation during a period of CAMo (Fig. j^b)). If, 
however, a cell is located in a central position, as in the 
case of 5 cells, this cell cannot establish a stable axis of 
internal polarization. It seems likely that this lack of ori¬ 
entation leads to the elevated instability we observed for 
CAMo states of such systems. 

To test these heuristic ideas we have developed a com¬ 
putational model [Si] generalizing the CPM [55[ [2S] to 
account for both internal cellular polarization and inter¬ 
cellular coupling. In the CPM, a cell is represented as a 
simply connected set of grid sites on a two-dimensional 
lattice, and thereby cell shape is explicitly represented. 
The model accounts for mechanical properties of cells and 
cell-cell adhesion. Previous generalizations of the CPM 
have implemented cell polarity and ensuing cell migration 
in a global fashion EiiEa upon adapting ideas from flock¬ 
ing models [18] : the overall polarity of a cell is described 
by a polarity vector, and it is assumed that there is a pos¬ 
itive feedback between a cell’s displacement and polarity. 
While these assumptions provide a simple and efficient 
way to model interactions between a cell and its mechan¬ 
ical environment, they do not resolve internal polariza¬ 
tion mechanisms. In fact, there are complex biochemical 
networks, including Rho family GTPases and membrane 
lipids, that regulate the assembly of the actin cytoskele- 
ton and thereby the formation of cell protrusions. Re¬ 
cently, computational models have been developed which 
couple rather sophisticated reaction-diffusion networks to 
the dynamics of membrane protrusions [ 321133 ] . These 
studies have provided important insights into the spa¬ 
tially resolved signaling processes within cells and how 




















































FIG. 4. (color), (a) Heat map of the relative nuclei positions with respect to a reference nucleus located at the lower border 
of the system from experimental and theoretical data (for details on the plot generation see section S8 [36]). (b) Schematic of 
possible polarization alignments during CAMo for different cell numbers, (c) For the CPM, mean magnitude of polarization 
and its error is plotted against the radial cell position r normalized by the maximal radial cell position ro. 


they are affected by cell shape. Here, in order to describe 
the dynamics off small cell groups, we used an intermedi¬ 
ate approach between flocking-type CPM models for cell 
assemblies EH 132] and detailed reaction-diffusion models 
for individual cells EH |33| . Specifically, in our compu¬ 
tation model we employed an internal polarization field 
within each individual cell to achieve the spatial resolu¬ 
tion of microscopic models. At the same time, the numer¬ 
ical algorithm is entirely rule-based (rather than based on 
complex reaction-diffusion networks) to retain the com¬ 
putational efficiency of CPMs. Furthermore, to account 
for the effects of cell-cell communication via mechano- 
transduction, the local dynamics of the internal polar¬ 
ization field is coupled to a cell’s membrane protrusions 
over a finite signaling range. This creates a positive feed¬ 
back loop integrating intracellular fluctuations and ex¬ 
ternal (mechanical) stimuli and gives rise to spontaneous 
cell polarization. To match the rotation statistics to the 
experiments, we simulated cells of fixed (average) size on 
circular islands at fixed cell density. We then performed a 
parameter sampling by varying cell adhesion, the range of 
intracellular mechanical signaling, and the strength of cy- 
toskeletal forces relative to contractile forces. For a more 
detailed and technical description of the model please re¬ 
fer to SI in the supplementary material |36j . which also 
contains a list of the model parameters used. 

The model reproduces the symmetry breaking into ro¬ 
tational states found by experiment (see |3S| for movies). 
Analyzing the numerically generated cell tracks analo¬ 
gously to experimental data, we found CAMo as well as 
DisMo (see section S7 [36]). Monte Carlo Steps were ad¬ 
justed to real time by matching the CAMo peak positions 
Ct]\f. We found the same steady decrease of Hjv with in¬ 
creasing cell number as in the experiments (Fig. ic)). 
Moreover, simulation data also exhibit an increase in 


CAMo persistence with increasing cell number for 2-,3- 
and 4-cell systems, and reproduce the discontinuity be¬ 
tween 4- and 5-cell systems (Fig. |^b)). (This feature is 
also observed when alternative measures for persistence 
are used; see S6 EH-) The discontinuity in persistence is 
accompanied by the same topological transition in cell ar¬ 
rangement as found by experiment (Fig.|^a) lower part). 
Assessing the characteristics of internal cell polarization 
in the model, we found a systematic decrease of the mean 
magnitude of the front-rear polarization with decreasing 
distance to the system center (Fig.|^c)). These findings 
clearly show that a cell in the center of the pattern is un¬ 
able to establish a stable axis of polarization and hence 
destabilizes collective behavior throughout the system. 

Here, we have presented a mesoscopic experimental 
setup in which the emergence and persistence of col¬ 
lective behavior is analyzed for a small number of cells 
in confined geometry. We showed that it is possible to 
obtain controlled migrational cell states, which may be 
classified as disordered and coherent angular motion. 
Both experiments and simulations showed consistently 
that persistence of the coherent state increases with 
the number of confined cells for small cell numbers but 
then drops abruptly in a system containing 5 cells. This 
is attributed to a geometric rearrangement of cells to 
a configuration with a central only weakly polarized 
cell. It reveals the decisive role of the interplay between 
local arrangement of neighboring cells and the internal 
cell polarization in collective migration. Future studies 
combining well-controlled cell assemblies confined to 
micropatterns and computational models may help to 
evaluate and characterize migrational phenotypes and 
identify mechanisms that play a key role in cell-to-cell 
mechano-transduction and finally in the emergence of 
collective behavior. 
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